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magnetic properties 
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The structural and magnetic properties of Fe octaethylporphyrin (OEP) molecules on Cu(OOl) 
have been investigated by means of density functional theory (DFT) methods and X-ray absorption 
spectroscopy The molecules have been adsorbed on the bare metal surface and on an oxygen- covered 
surface, which shows a \/2 x 2v^i?45° reconstruction. In order to allow for a direct comparison 
between magnetic moments obtained from sum-rule analysis and DFT we calculate the dipolar term 
7(Tz), which is also important in view of the magnetic anisotropy of the molecule. The measured X- 
ray magnetic circular dichroism shows a strong dependence on the photon incidence angle, which we 
could relate to a huge value of 7{Tz), e.g. on Cu(OOl) 7{Tz) amounts to -2.07 /xb for normal incidence 
leading to a reduction of the effective spin moment ms-\-7{Tz). Calculations have also been performed 
to study the influence of possible ligands such as CI and O atoms on the magnetic properties of 
the molecule and the interaction between molecule and surface, because the experimental spectra 
display a clear dependence on the ligand, which is used to stabilize the molecule in the gas phase. 
Both types of ligands weaken the hybridization between surface and porphyrin molecule and change 
the magnetic spin state of the molecule, but the changes in the X-ray absorption are clearly related 
to residual CI ligands. 



I. INTRODUCTION 

Hybrid systems of paramagnetic molecules and metal- 
lic surfaces are new promising materials for magneto- 
electronic devices and may in future replace devices 
with artificial nano- structures.^ ^ The ongoing search for 
new materials which can be used in magnetic switch- 
ing devices has brought phthalocyanine and porphyrin 
molecules with transition metal center into focus, be- 
cause they arrange themselves flat on surfaces and are 
relatively easy to handle experimentally under UHV con- 
ditions. Recent experiments with Co (II) tetraphenylpor- 
phyrin (CoTPP) molecules on Ni/Cu(001) have shown 
that the spin state of the Co atom can be controlled 
by NO ligands, which are attached or removed by tem- 
perature treatment.^ Miguel et al. tailored the coupling 
between Fe(II) octaethylporphyrin (Fe OEP) molecules 
and an oxygen covered Ni/Cu(001) substrate by the same 
type of ligands.^ It has been observed that in the pres- 
ence of an intermediate oxygen layer the coupling be- 
tween the paramagnetic molecules and the ferromagnetic 
(FM) surface is of antiferromagnetic nature, whereas the 
investigated systems couple ferromagnetically without 
oxygen Theoretical investigations have revealed that 
Fe OEP molecules are chemisorbed at ferromagnetic sur- 
faces whereby the magnetic coupling is mediated by the 
N-atoms.^ 

Transition metal-based porphyrin molecules have been 




FIG. 1: Schematic structure of the Fe OEP molecule. Hy- 
drogen atoms are not shown. For the calculations the ethyl 
groups CHs are replaced by H atoms. Ligands are denoted 
by L and stick out of the molecular plane. 



studied intensively on ferromagnetic substrates fo- 
cussing on the magnetic coupling between surface and 
molecule'^HI^ and the influence of the anisotropy of the 
substrate on the magnetic properties of the molecule. 
Nonmagnetic substrates have been used to study spin and 
orbital magnetic anisotropiesl^^^ as well as correlation 
and hybridization effects with the substrate. Here, 
we use a non-magnetic substrate to unravel the magnetic 
anisotropy of the molecule and it's dependence on ligands 
and surface oxidation, whereby special focus is on the im- 
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FIG. 2: (Color online) (a) \f2 x 2v/2R45° missing-row reconstruction of the Cu(OOl) surface with 0.5 ML O. Cu atoms are 
shown in different gray shades (light gray corresponds to the surface layer and black to the fixed bottom layer.) The smaller 
(red) circles indicate the O atoms. Numbers mark the different positions of the Fe center of the FeP, namely missing-row (1), 
on top of O (2), and hollow site (3) position, (b) Charge density plot of the reconstructed O/Cu(001) surface in comparison 
with the plain Cu(OOl) surface (c). Dark (red) circles correspond to the charge distribution of the Cu atoms, lighter (orange) 
spots in (b) indicate the charge of the surface 0-atoms. The bright (greenish) areas in (b) (marked by the dotted line) show 
the charge depletion due to the missing-row reconstruction. 



portance of the dipole term l{Tz). We address the issue 
of how the magnetic properties like magnetic anisotropy, 
and the hybridization between molecule and substrate 
change if the FM film is absent, i.e., we study the ad- 
sorption of Fe porphyrin molecule on a bare Cu(OOl) sub- 
strate and pursue the question whether an 0-interlayer 
has an influence on the magnetic properties as in the pres- 
ence of an FM film. This has been done by performing 
X-ray absorption measurements and density functional 
theory (DFT) calculations using the VASP code.^^ For a 
proper description of molecular bonding van der Waals 
forces are important. However, they are not included in 
the standard DFT description, therefore, we have used 
the semi-empirical form of Grimme to account for van 
der Waals interaction. The use of a fully ah initio van 
der Waals term is not feasible for large systems as in the 
present case. 

Parallel to the present ah initio investigations X-ray 
absorption experiments have been carried out for Fe OEP 
on plain Cu(OOl) and in presence of an oxygen interlayer, 
which goes along with a x 2v/2R45° reconstruction of 
the Cu surface, see Fig.[2](a). We have performed angle- 
dependent X-ray absorption spectroscopy (XAS) and X- 
ray magnetic circular dichroism (XMCD) experiments in 
high magnetic fields at low temperatures. From these 
measurements we can evaluate the effective spin moments 
via sum-rules derived by Carra et al.l^ In order to com- 
pare the effective spin moments extracted from experi- 
ment and the calculated spin moments we determine the 
dipolar term which connects the two quantities according 
to the work of van der Laan.l^ 

In the present DFT calculations the Fe OEP molecules 
are modeled without the outer ethyl groups and are passi- 
vated with H atoms only, referenced as FeP in the follow- 



ing, see Fig.[T] Furthermore, in the experiment Fe OEP 
molecules are usually decorated with additional ligands 
to stabilize them. In our case we use two different types 
of ligands: pyridine (C5H5N named Py in the follow- 
ing) , which decouples from the molecule during degassing 
of the molecular powder in vacuum prior to the sample 
preparation which yields Fe(II)OEP in its pure form and 
atomic CI for which scanning tunneling microscopy mea- 
surements indicate that the CI may partially remain in 
the system. We will show that the choice of the ligand, 
which stabilizes the molecule in air has a drastic influence 
on the magnetic properties of the molecule. 

The paper is organized as follows. After a brief presen- 
tation of the theoretical methods including a discussion 
of the preparation of the surfaces in Secjllj in Sec. [Hi 
the experimental set up and sample preparation will be 
explained. A discussion of the results can be found in 
Secs.[TV| and [V| where magnetic and structural properties 
of the FeP on the two different surface structures — plain 
Cu(OOl) and oxidized Cu(OOl) surface — as well as the 
influence of CI and O ligands will be considered. Finally, 
a conclusion and an outlook are given in the last section 
(Sec.lvnl). 



II. COMPUTATIONAL METHODS AND 
SURFACE RECONSTRUCTION 

We have performed density functional theory (DFT) 
calculations employing the VASP code^ with the pro- 
jected augmented wave potential (PAW) method. The 
exchange correlation functional has been desc ribed via 
the generalized gradient approximation (PW91) plus 
a Hubbard- /7 correction in order to incorporate the lo- 
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calized character of the d orbitals of the Fe-atom in 
the molecule and to get a reasonable HOMO-LUMO 
gap. Here, we make use of the method of Dudarevl^, 
where the Coulomb (/7) and exchange (J) interaction 
enter in the Hamiltonian as an effective value, whereby 
U — J = 3.0 eV has been used for the Fe d orbitals. 
The value has been previously shown to describe FeP 
properly. However, the description of molecular sys- 
tems in DFT does not only suffer from the failure of 
the description of correlation effects, but also does not 
account for van der Waals forces. Recently, several at- 
tempts have been made to include van der Waals in- 
teraction in DFT.^^ For such large systems as in 
the present study we use the semi-empirical method of 
Grimme as implemented in VASP. Fully ab initio meth- 
ods are to expensive for the systems under consideration. 
However, a semi-empirical approach should already sig- 
nificantly improve the description of the molecular bond- 
ing. All calculations have been performed for the F-point 
and a cut-off energy of 400 eV for the plane- waves. In or- 
der to study the magnetic and electronic properties of 
FeP on Cu(OOl) with and without an intermediate oxy- 
gen layer we have used an 8 x 8 lateral unit cell and 3 
monolayer (ML) of Cu(OOl). The bottom layer is kept 
fixed to simulate bulk-like behavior. The rest of the sys- 
tem is fully relaxed and the vacuum size is about 14 A. 

If Cu(OOl) is covered by 0.5 ML of oxygen this leads 
to a V2 X 2V^R45° surface reconstruction, see Fig.|| 
(a).^^^^ After relaxation the O atoms are embedded al- 
most in the surface layer being only 0.16 A above the 
Cu atoms. This is in good agreement with results from 
low energy electron diffraction (LEED) measurements 
by Zeng et al?^ This reconstruction has significant in- 
fluence on the electronic structure of the surface layer. 
Figure[2] shows the charge distribution of the Cu(OOl) 
X 2V^R45° surface (b) in comparison to the plain 
Cu surface (c). The charge depletion along the missing- 
rows (blue-green areas) is clearly visible. Due to the 
missing-row reconstruction the electronic structure of the 
0-terminated surface is more complex compared to the 
plain Cu(OOl) surface and offers a number of additional 
adsorption positions for the molecule, e.g., the Fe sit- 
ting on the missing-row or on top of an oxygen atom, cf. 
Fig-ia). 

The surface or the ligands may not only influence 
the magnetic spin state but also the magneto-crystalline 
anisotropy (MCA), i.e., the energy needed to switch the 
magnetic orientation from the easy axis to the hard axis. 
Here, we will consider only the spin-orbit coupling in- 
duced magnetic anisotropy, for details see Ref. 29, The 
perturbation of the original system-originating spin-orbit 
coupling can be written as 

Hso = C(OL • S = ar){LA + LySy + L,S,) (1) 

4c^r or 

with C being the spin-orbit coupling strength. For 3(i 
transition metals, ^ is in the range of 50-70 meV which is 



small compared to the usual band width and hence, the 
approximation is reasonable. Since the first order term 
vanishes, the first significant contribution comes from the 
2nd order corrections. For 3(i metals, higher order terms 
are ignored. Therefore, the spin-orbit contribution to the 
energy, from the 2nd order perturbation theory reads 

^SO = -C yriuTlo . [6) 

Here \o) and \u) correspond to occupied and unoccupied 
states, respectively, in the basis of |/m,<j), where l^m^cr 
are orbital, magnetic, and spin quantum numbers, re- 
spectively. Eu, riu and Eo, Uq are energy eigenvalues and 
occupations of unoccupied and occupied levels. The elec- 
tronic states of our interest are antibonding states of Fe-N 
p — d hybridization arising mostly from Fe d states in FeP. 
The contribution to the MCA is inversely proportional to 
AE = Eu — Eo (see Eq.[3|. Thus states close to the Fermi 
energy (^f) are the most important ones. The eigenval- 
ues are taken from our calculated orbital-projected den- 
sity of states and the 2nd order perturbation has been 
calculated for magnetization axis along (001) and (100), 
i.e out-of-plane and in-plane directions. 

III. EXPERIMENTAL DETAILS 

Iron OEP was synthesized by reduction of Fe OEP(Cl) 
using hydrazine and isolated as its dipyridine complex Fe 
0EP(Py)2.^ Differential Scanning Calorimetry (DSC) 
studies of the isolated material showed loss of the pyridine 
ligands at 146-154 ^C.l^ No further transitions were de- 
tected up to 350 °C. Paramagnetic Fe OEP was deposited 
in UHV on Cu(OOl) and on an oxygen-reconstructed 
Cu(OOl) surface (held at room temperature) from a tan- 
talum crucible at about 200 °C. Prior to the deposition, 
Fe 0EP(Py)2 was degased in UHV at 200 °C for about 
5 h to ensure the loss of the pyridine ligands. After depo- 
sition, N K edge XA spectra displayed no tt* resonance 
in the energy range from 397 eV to 403 eV in addition 




FIG. 3: (Color online) Sketch of the experimental setup. The 
molecules (ovals) are placed on the Cu(OOl) surface. The 
X-ray absorption measurements are carried out for different 
angles between the surface normal and the photon beam 
(k). The magnetic field is always oriented parallel to the k- 
vector of the photon beam and the electric field, which probes 
the d orbitals, is in the plane of incidence, perpendicular to 
the photon beam (p-polarized). 



4 





to the resonances of the porphyrin macrocycle, showing 
that pyridine is not present on the surface. For the 2nd 
set of samples Fe OEP(Cl) was purchased from Porphyrin 
Systems and deposited from a tantalum crucible at about 
200 °C without further treatment. Angle dependent XAS 
and XMCD measurements have been performed at the 
beam lines UE46 PGMl at BESSY II (samples with Py) 
and ID08 at the ESRF (samples with CI) in an applied 
magnetic field (temperature) of 5.9 T (5K) at BESSY II 
and 5.0 T (8K) ESRF, respectively. The XAS signal was 
detected by means of total electron yield and normalized 
to a reference signal upstream to the experiment. The k 
vector of the X rays was aligned parallel to the magnetic 
field for different incidence angles 6>, see Fig. [3] The third 
(first) harmonic of the undulator with an energy resolu- 
tion set to 150 meV (250 meV) and a circular polarization 
degree of 85% (100%) was used for the measurements per- 
formed at BESSY II (ESRF). Coverages of 0.4 ML were 
chosen to ensure that only molecules that are in direct 
contact with the substrates contribute to the XAS sig- 
nal. In Fig.|4 scanning tunneling microscopy (STM) im- 
ages of 0.4 ML of Fe OEP(Cl) on Cu(OOl) taken directly 
before the XAS measurements at ID08 are shown. It is 
obvious that the porphyrin molecules adsorb fiat on the 
surface. Bright and dark centers of the molecules may 
be attributed to CI ligands on top of the iron centers 
and to the Fe OEP where the CI ligand has desorbed, 
respectively. 




IV. EXPERIMENTAL RESULTS 

In Fig.[5]the integrated Fe L3 XMCD signal of 0.4 ML 
Fe OEP on Cu(OOl) and oxygen- covered Cu(OOl) is plot- 
ted as a function of applied magnetic field along 0° and 
55° at T = 5K. In both cases the XMCD at 55° is sig- 
nificantly higher than the one at 0° . This can be a result 
of either an in-plane magnetic anisotropy, a large mag- 
netic dipolar term (T^), or a mixture of both effects. A 
detailed discussion is found in Sec. lVI A"2l For Fe OEP 
on Cu(OOl) (Figjsja)), the magnetization curves are al- 
most straight indicating that the Fe ions are magneti- 
cally not saturated at 5 = 6 T. The bended curves for Fe 
OEP on O/Cu(001) (Fig. |5|b)) on the other hand sug- 
gest that magnetic saturation can be reached at 55° and 
5 = 6T for this system. The Fe ^2,3 XAS and XMCD 
signals of 0.4 ML Fe OEP (Fe OEP (CI)) on Cu(OOl) 
and oxygen- covered Cu(OOl) are shown in Figs.[6]^a) to 
^d) Qe) to|6jf)) for three different incidence angles at 
B = 5.9 T {B = 5.0 T) and a temperature of T = 5K 
(T = 8K). The Fe L3 XAS signal of Fe OEP on Cu(OOl) 
comprises of a main resonance at 708 eV mostly visible 
for normal incidence and three shoulder-like features at 
705.7 eV, 706.6 eV, and 707.2 eV. The integrated Fe ^2,3 
XAS signal displays a strong angle-dependent variation 
even for circularly polarized light. For linearly polarized 
X rays, a ratio of 1.5 between the integrated Fe ^2,3 XAS 
signal (not shown) for 0° and 65° incidence angle is found 



FIG. 4: (Color online) STM images of 0.4 ML of Fe OEP(Cl) 
on Cu(lOO) at room temperature (a): 1.3V, 0.08 nA; (b): 
1.2V, 0.17nA. 



indicative for an Fe 3d hole density distribution that is 
predominantly oriented in the surface plane. The Fe I/3 
XAS signal of Fe OEP on O/Cu(001) displays two main 
resonances at about 707.5 eV and 709 eV. The overall Fe 
L3 XAS intensity is shifted by approximately 1.0 eV to 
higher photon energy compared to Fe OEP on Cu(OOl). 
This shift can be attributed to a charge transfer from 
the Fe atom in the molecule to the oxygen atom in the 
substrate. This is an indication that the valence of the 
Fe atom is trivalent on the oxygen covered surface. Sim- 
ilar trends have been detected for oxygen- covered mag- 
netic surfaces.l^ The integrated Fe ^2,3 XAS signal dis- 
plays little variation with the incidence angle represent- 
ing a nearly isotropic 3d hole density distribution. The Fe 
L2,3 XMCD signal of Fe OEP on Cu(OOl) is smaller than 
the one on O/Cu(001). The Fe L3 XMCD signal of Fe 
OEP on Cu(OOl) (O/Cu(001)) displays a double peak at 
706.8 eV and 707.8 eV (707.4 eV and 709 eV) for 55° and 
65° incidence angle and a resonance at 708.1 eV (709 eV) 
for 6 = 0° that is smaller by about a factor 4 for both sys- 
tems. Since except for the measurements of Fe OEP on 
O/Cu(001) at ^ = 55° and 65° the Fe magnetic moments 
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were not fully aligned by the magnetic field {cf. Fig.[5| 
the strong angle-dependent variation of the XMCD signal 
can be a result of either an in-plane magnetic anisotropy 
or a large magnetic dipolar term (T^), i.e. , an ani sotropy 
of the spin-density distribution {cf. Sec. VI A 2). As it 
will be discussed in section [VI A 2| the prominent angular 
dependence is due to a strong angular-dependent con- 
tribution and the magnetic anisotropy of the molecules 
is very small. 

The Fe ^2,3 XAS and XMCD signals of 0.4 ML Fe 
OEP(Cl) on Cu(OOl) (Fig.jefe) and (f)) display an over- 
ah similarity to the ones of 0.4 ML Fe OEP on Cu(OOl) 
(Fig.|6|a) and (b)) for the spectra measured at grazing 
incidence and the magic angle with a less pronounced 
splitting of the double peak feature of the Fe L3 XMCD. 
On the other hand the XAS and XMCD signal for normal 
incidence are markedly different. Most prominently the 
Fe L2,3 XMCD of 0.4 ML Fe OEP(Cl) on Cu(OOl) is big- 
ger by a factor of about 4 than the one of 0.4 ML Fe OEP 
on Cu(OOl) at normal incidence resulting in an integrated 
XMCD signal that is basically the same for the different 
measurement angles. This shows that the electronic and 
magnetic properties of Fe OEP(Cl) on Cu(OOl) are differ- 
ent from the ones of 0.4 ML Fe OEP on Cu(OOl). Since 
the STM pictures presented in Fig. [4] show bright and 
dark molecular centers with a ratio of about 1:3 it can 
be assumed that part of the molecules have a CI ion at- 
tached to their iron center and others not. In this case 
the spectra presented in Fig.|6|e) and (f) are superposi- 
tions of two contributions that cannot be disentangled 
easily. The Fe ^2,3 XAS and XMCD signals of 0.4 ML Fe 
OEP(Cl) on O/Cu(001) (Fig.jefg) and (h)) present ba- 
sically the same spectral features and angle dependence 
than the ones of 0.4 ML Fe OEP on O/Cu(001) (Fig.[6];c) 
and (d)), showing that the electronic and magnetic prop- 
erties are comparable even though CI ions may be present 
on the surface for the former sample. The effect of CI 
attached to the metal center of the molecules will be dis- 
cussed in SectionlVll The results of the sum rule anal- 
ysis of the spectra shown in Fig. [6] is presented later on 
in this manuscript together with the theoretical results. 
The number of d holes has been assumed to be n/^ = 4 
for the molecules on the plain Cu(OOl) surface, which is 
in good agreement with the theoretical values obtained 
from the present DFT calculations {uh = 3.95). In case 
of the oxidized surface n/^ = 5 has been used, because 
the ^2,3 edge spectra of Fe are shifted by about 1 eV to 
higher energies, which can be assessed as an indication 
for Fe^+.EI 



V. FeP ON O/Cu(001) AND Cu(OOl) 

A. Calculated structural properties 

Supercell calculations have been performed to deter- 
mine the ground state geometry of the Fe porphyrin 
molecule on the plain metal and the oxidized surface. 
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FIG. 5: (Color online) Integrated Fe L3 XMCD signals of 
0.4 ML of Fe OEP on Cu(OOl) (a) and on O/Cu(001) (b) at 
T = 5 K as a function of applied magnetic field aligned parallel 
with the k vector of the X rays along ^ = 0° and = 55° . 



where different adsorption sites as well as different ori- 
entation relative to the substrate have been taken into 
account. The results are summarised in Tab. HI 



1. FeP on Cu(OOl) 

In case of a clean Cu(OOl) substrate the FeP molecule 
prefers the hollow site position, whereby the four N atoms 
are located on top of the neighboring Cu atoms of the 
surface layer, see Fig.[7| (a). The distance between the 
central Fe atom and the surface amounts to 3.10 A, if no 
vdW interaction is considered. Taking into account van 
der Waals forces in the approximation of Grimme the 
distance between molecule and surface is significantly re- 
duced by « 0.4 A to 2.66 A. This bonding distance on 
the nonmagnetic substrate is even with van der Waals 
interaction much larger compared to bonding distances 
obtained for such molecules on FM surfaces i.e., FeP 
or Fe phthalocyanine molecules are chemisorbed on mag- 
netic surfaces. This means that in contrast to a magnetic 
metal surface the hybridisation between molecule and the 
bare Cu(OOl) surface is much weaker. A rotation of the 
porphyrin molecule by 45° such that the N atoms are 
now located between the Cu atoms of the surface layer 
is connected with an increase of the energy by 0.311 eV 
(without van der Waals 0.068 eV). The predilection of 
the N atoms for the Cu atoms of the surface is related 
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FIG. 6: (Color online) Fe L2,3 XAS ((a), (c), (e), (g)) and XMCD ((b), (d), (f), (h)) signals of 0.4 ML Fe OEP on Cu(OOl) 
and (O)/Cu(001) substrates (panel (a) to (d)) and Fe OEP CI (panel (e) to (h)) on substrates of Cu(OOl) and (O)/Cu(001) 
for different photon incidence angle ^ at T = 5K ((a)-(d)) and T =8K ((e)-(h)), respectively. Note the different scales for the 
XMCD signals. 



to the fact that the hybridization between molecule and 
surface is mediated by the N-atoms, cf. Sec. |VB 3) This 
effect has been previously observed for porphyrin and ph- 
thalocyanine molecules on ferromagnetic transition metal 
surfaces!^ All other positions with Fe on top of a Cu 
atom or on a bridge position turned out to have signif- 
icantly higher energies, see Tab.[Tj In case of the bridge 
position the energy difference to the hollow site position 
amounts up to 5 eV even without van der Waals interac- 
tion. For Fe on top of Cu the energy is approximately 
0.5 eV higher than for the ground state, independent of 
whether the N-atoms are on top of Cu or not. These 



configurations will not be considered further. 



2. FeP on O/Cu(001) 

The deposition of oxygen on Cu(OOl) is accompanied 
by a y/2 x 2^/2 R45° reconstruction of the surface and 
therefore, the surface offers adsorption sites being differ- 
ent from the ones on Cu(OOl), see Fig.[2|a). However, 
due to the missing-rows only few positions are available 
which fit to the 4-fold symmetry of the molecule. Since, 
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TABLE I: Calculated energy differences AE and Fe-surface distances dpe-s for FeP on different adsorption positions of bare 
Cu(OOl) and >/2 x 2a/2R45° O/Cu(001), with respect to the total energy of the ground state configuration. The values obtained 
without van der Waals correction are given in brackets. 
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FIG. 7: (Color online) Schematic view of the ground-state 
configurations of FeP on Cu(OOl) (a) and in presence of an 
intermediate oxygen layer and a V2x 2v^ R45° reconstructed 
surface (b). Gray spheres mark the Cu atoms of the substrate, 
purple corresponds to O atoms. The FeP consists of an Fe 
central atom (red) surrounded by 4 N-atoms (green) and 20 
C atoms (yellow). The ethyl groups have been replaced by 12 
H atoms (cyan). 



O is attractive for Fe in addition to the previously inves- 
tigated adsorption sites, configurations with Fe on top 
of O and on the missing-row have been taken into ac- 
count. From the calculations it turned out that the most 
stable structure is indeed the missing-row position, i.e., 
the Fe atom sits directly over the missing-row with the 
4 N atoms on top of Cu atoms, see Fig.[7| (b). Here, the 
Fe atom has two O neighbors in (-110) direction. This 
leads to a tiny deformation of the N cage, i.e. the an- 
gle between Fe and the two N atoms with O amounts to 
89.4°. The hollow site position, which was the favorable 
position in case of the plain Cu(OOl) surface, is without 
van der Waals interaction only 15 meV higher in energy. 
Configurations with Fe on top of an O atom or on a 
bridge-like position (Fe in the middle between Cu and 
O) can be found about 40meV (without van der Waals 
forces) above the ground state configuration, whereby for 
Fe on top of O only three of the four N atoms have a Cu 
counterpart to hybridize with. In case of the less sym- 
metric bridge-like position none of the N-atoms has the 
chance to interact directly with a Cu atom, instead two 
of the C atoms of each pyrrole unit sit on top of the Cu 
atoms. However, the relatively small energy differences 
between the different positions suggest that all four po- 
sitions may be partially occupied at finite temperatures. 
The picture drastically changes if van der Waals forces 



are included. Although the ground state configuration 
remains the same and also the energetic sequence of the 
adsorption positions is the same, the missing-row position 
becomes more stable compared to the hollow site, top of 
oxygen, and bridge-like positions, see Tab.jlj The missing- 
row position is now 0.566 eV lower in energy compared to 
the hollow site position. Obviously, the long-range van 
der Waals forces drift apart the energies of the differ- 
ent adsorption sites, which is related to the reduction of 
the molecule-surface distance and the related stronger 
surface-molecule interaction. The results for the 45°- 
rotated positions are not given in Tab.jlj because they 
are the less favorable configuration on each adsorption 
site. For example rotating the molecule on the missing- 
row position by 45°, such that two of the N-atoms end 
up sitting on top of O, costs approximately 0.7 eV with- 
out van der Waals correction, and the magnetic moment 
vanishes. Without van der Waals forces the distance be- 
tween molecule and surface is so large (> 3.2 A) that the 
interaction with the surface is rather weak, whereas with 
van der Waals forces the hybridization between molecule 
and surface increases. The distance between Fe atom 
and oxygen is now reduced to 2.65 A for the ground state 
configuration. As a consequence, asymmetric configura- 
tions will be rather unlikely, e.g., the configuration with 
Fe on top of O is now 0.580 eV higher in energy because 
only three N atoms have the possibility to hybridize with 
Cu atoms of the underlying surface, whereas the fourth 
N-atom is located on the missing-row position, see Fig.[7| 
(b). This clearly demonstrates the importance of van der 
Waals corrections in this type of systems. However, in 
experiment other than missing-row position may be par- 
tially occupied for larger molecule concentrations and in 
case of low growth temperatures, where the mobility of 
the molecules is limited and they may stick to unfavor- 
able positions. Therefore, in the following we discuss the 
ground state configuration and the configuration being 
closest in energy to the ground state. It should be noted 
that the inclusion of van der Waals forces is responsible 
for a roughening of the surface layer. The surface, which 
had previously been mainly flat, shows now a distinct 
buckling, i.e., the height of the O atoms varies by 0.78 A. 
O atoms under the molecule move inwards whereas the 
other stick out of the surface. The same effect is ob- 
served on the plain Cu(OOl) surface, but it is much less 
expressed. 
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FIG. 8: (Color online) Calculated density of states of Fe d orbitals in FeP/Cu(001). Arrows denote the spin-up and spin-down 
configuration, respectively. The Fe atom is adsorbed on the hollow site position with N on top of Cu (a) and (c). In case 
of (b) and (d) the molecule is rotated by 45°. The results obtained from calculations with (without) van der Waals forces in 
the approximation of Grimme are given in the top row (bottom row). Differently colored lines denote different d orbitals as 
indicated in the legend in (a). 



B. Calculated magnetic properties 

From here on we focus only on the two most stable 
configurations, i.e., the ground state and state which is 
closest in energy to the ground state. 



1. FeP on Cu(OOl) 

The FeP molecule on Cu(OOl) carries a spin moment 
of 2 /iB and is in an intermediate spin state 5* = 1 for all 



possible adsorption sites, which agrees with our results 
from the sum-rule analysis. For a more direct compari- 
son of the magnetic moments with the effective magnetic 
moments from the sum-rule analysis the dipolar term has 
to be included, see discussion in Sec. |VI The = 1 
spin-state turns out to be very stable and is not changed 
when van der Waals forces are included. If in addition 
a high-spin (4/iB) state can be stabilized it has always 
much higher energy, e.g. for the ground state configu- 
ration the high-spin state lies without van der Waals in- 
teraction ?^0.5eV above the intermediate spin-state. The 
density of states of the Fe d orbitals of FeP on the hollow 
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FIG. 9: (Color online) Calculated density of states of Fe in FeP/ s/2 x 2^R45°O/Cu(001). Arrows denote the spin-up and 
spin-down configuration, respectively. The Fe atom is placed on the hollow site ((a), (c)) and the missing-row position ((b) and 
(d)) with the N-atoms on top of Cu. Results including van der Waals (vdW) forces are shown in the top row, results without 
van der Waals interaction are given in the bottom row. Differently colored lines represent different d orbitals as indicated in 
the legend in (a). 



site position can be found in Fig.jS] For the ground state 
hollow site position (N on top of Cu) and the 45° rotated 
molecule (N between Cu) a (dxy)^ ^ (<^7r)^, {^'^z'^Y c^^" 
figuration appears if van der Waals forces are included. 
This corresponds to recent findings for the gas phase 
FeP molecule In contrast to that, without van der 
Waals interaction the d^i orbital is mostly filled, such 
that the occupancy corresponds to (d^y)^ ^ (<^7r)^, {dz^)'^^ 
see Fig. [8) Besides the change in occupation of the z- 
oriented orbitals the density of states with and without 
van der Waals interaction are very similar. In contrast 
to the gas phase FeP molecule the density of states is 
broadened and no clear HOMO-LUMO gap can be ob- 



served. The ligand field splitting which can be obtained 
from the energy difference between the dxy and dx2_y2 
orbital amounts to ~ 1 eV in all cases, when the artificial 
splitting due to the Hubbard-like U — J term has been 
subtracted. 



2. FeP on O/Cu(001) 



As discussed in Sec. |V A l] the adsorption site of the 
molecule is different on vTx 2V2R45°O/Cu(001) com- 
pared to Cu(OOl), however this has just as little influ- 
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ence on the magnetic structure resulting in a slightly 
deformed Fe-N cage. The magnetic configuration re- 
mains unchanged, i.e., the Fe atom is in an intermedi- 
ate spin state and carries a spin moment of 2 jub- For 
the Fe atom on the missing-row position, the occupa- 
tion of the orbitals is the same as for the previous case 
{{dxy)'^^ (<^7r)^7 {dz^)^)^ independent from whether van 
der Waals interaction is included or not, see Fig.[9]^b), 
(d). The HOMO-LUMO gap amounts to 0.9 eV, which 
is only slightly smaller than the value obtained for the 
gas phase molecule^^. In contrast to the observations on 
the missing-row position, for molecules on the hollow site 
position the orbital is completely filled, see Fig.|9ja) 
and (c). Interestingly, the system has - at least in the 
unoccupied part - narrow energy bands, i.e., sharp peaks 
in the density of states, which is different from the ob- 
servations on the plain Cu(OOl) surface. Here, the shape 
of the density of states is more alike to the gas phase 
molecule Only the occupied d-j^ and d^^ orbitals are 
somewhat broadened compared to the free Fe porphyrin, 
which stems from interaction with the surface, see dis- 
cussion in the next section. This broadening is not ob- 
served for the hollow site position without van der Waals 
interaction (Fig.[9]^c)), which is explained by the huge 
molecule-surface distance of 3.48 A. 
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3. Hybridization and coupling to the surface layer 

So far the electronic, magnetic and structural proper- 
ties of the molecules themselves have been studied, here 
we focus on the interaction between molecule and surface. 
It has been shown in Sec. lV Al that the distance between 
molecule and surface even with van der Waals correction 
is about 2.7 A, thus the hybridization effects with the 
surface are weaker compared to adsorption on magnetic 
substrates.*^ Without oxygen the interaction between Fe 
and Cu atoms is small. On one hand small hybridiza- 
tion can be assumed from peaks around 2-3 eV below the 
Fermi level in the majority-spin channel, which appear 
for both systems, but no significant hybridization could 
be observed for the minority-spins, see Fig. 10 ^a). This 



suggests that the direct interaction between Cu atoms 
and Fe atom plays only a minor role. On the other hand 
as can be seen from Fig.p^a), without oxygen interlayer 
the positions of the N orbitals overlap with Cu and Fe or- 
bitals, i.e., the interaction between molecule and surface 
is also mediated by the N atoms. This matches with the 
fact that the preferred orientation of the molecule on the 
surface is with the N-atoms on top of Cu, cf. Sec. |V A| 
This type of indirect coupling has been observed on mag- 
netic surfaces, e.g. for FeP on Co(OOl).^ The charge of 
the Fe atom and the four N atoms is reduced by 0.0 le 
compared to the gas phase and the total charge of the Cu 
atoms in the vicinity of the molecule is increased com- 
pared to the free Cu(OOl) surface. 

In case of the oxygen-covered surface the situation be- 
comes more complex, see Fig.[lQ[b). Here, the hybridiza- 



FIG. 10: (Color online) Calculated density of states of Cu and 
O from the surface layer and density of states of N atoms and 
Fe d orbitals from the FeP molecule, (a) FeP on the hollow site 
position on the plain Cu(OOl) surface and (b) on missing-row 
position of the oxidized surface, respectively. Calculations are 
performed with van der Waals interaction.^^ Arrows denote 
the spin-up and spin-down configuration, respectively. 



tion between Cu and Fe d electrons seems to be somewhat 
larger compared to the plain Cu surface, especially in the 
majority channel the Cu d state positions coincide with 
the Fe d orbital. This means in contrast to the previous 
case the coupling is not dominated by the indirect cou- 
pling of Fe to Cu via N. Instead the Fe atom hybridizes 
with the neighboring O atoms, for example in the minor- 
ity channel at 1 eV below and 2 eV above the Fermi level. 
Similar as before the charge of the Fe atom is slightly re- 
duced compared to the charge of the Fe atom in the gas 
phase molecule. Here, a small charge transfer from the 
Fe atom to the O and Cu atoms under and next to the 
molecule is likely, because the oxygen charge is larger by 
0.014e as without the molecule. The same holds for the 
Cu charge. The N charge remains unchanged. Although, 
the changes of the individual charges are small, they are 
in line with the experimental findings for Fe OEP on the 
O covered surface, cf. SeclIVI However, the calculated 



charge variation is much too small to indicate the ex- 
istence of trivalent Fe. The fact that a charge transfer 
occurs from Fe to O indicates that the Fe atom directly 
hybridizes with the surface. Similar effects have been 
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FIG. 11: (Color online) Sketch of optimized structure of FeP 
with CI ligand on ^2 x 2^2 R45°O/Cu(001). 



observed for FeP on 0/00(001)1^. However, there the 
driving force is the magnetic exchange coupling which 
does not exist here. However, we observe a hybridization 
between O and N orbit als {s and p type) 4 eV and more 
below the Fermi energy, see Fig.[Toj which means that an 
indirect binding exists at least partially. 

Although there occurs a small charge transfer between 
molecule and surface and hybridization between the or- 
bitals of the molecule and the ones from the surface exist, 
the situation is completely different from the scenario on 
the magnetic substrates, where a magnetic coupling be- 
tween molecule and surface exists.— 



VI. INFLUENCE OF O AND CI LIGANDS 

Fe porphyrin molecules are delivered with stabilizing 
ligands such as 01 or pyridine (O5H5N). Although the 
ligands should be removed during the preparation pro- 
cess, especially 01 atoms may remain in the system and 
stick to the FeP molecules. This is suggested by STM 
and XAS measurements of Fe OEP(Ol) on Ou(OOl) and 
discussed in Sec. |IV[ In case of an oxygen interlayer there 
may also be some residual oxygen gas remaining in the 
chamber. Therefore, we have investigated the influence 
of these two types of ligands attached to the FeP molecule 
on the oxidized surface, respectively. In case of Ou(OOl) 
only 01 ligands have been considered, because no resid- 
ual oxygen gas is expected in this case. Single 01 or O 
atoms have been attached to the Fe-atom at the remote 
side from the surface, see Fig.[TT] Here, the discussion is 
restricted to calculations including van der Waals forces. 

Both types of ligands lead to a wide difference in elec- 
tronic and magnetic structure compared to the ideal FeP 
molecule on Ou(OOl) and O/Ou(001). Chlorine ligands 
provoke a destabilization of the molecule-surface hybrid 
system and lead to a concomitant increase of the dis- 
tance between surface and FeP molecule by 0.43 A for 
FeP/O/Ou(001), see TablejllJ No such destabilization 
could be observed on the plain Ou(OOl) surface, because 
it is related to the electrostatic repulsion between 01 and 
O atoms. On Ou(OOl) the distance between surface and 
molecule remains unchanged compared to the ligand-free 
case. Due to the relatively large distance to the surface. 



the hybridization with the surface atoms is comparable to 
the one without 01. In contrast to that, the distance be- 
tween the 01 and Fe atom amounts to 2.26 A (2.29 A on 
Ou(OOl)), which is of the same order as the bond lengths 
in the molecule itself. Consequently, the d states of Fe 
which stick out of the molecular plane (^^2, d^) hybridize 
with the 01 p states, and the unoccupied d states are 
shifted closer to the Fermi level, see FigjT2] Due to the 
weak coupling to the surface, the magnetic and electronic 
structure of the molecule is strongly influenced by the lig- 
and. If 01 is attached to the Fe atom, the spin moment 
increases from 1.99 to 2.62 /iB for the molecule on the ox- 
idized surface and from 1.93 to 2.50 /iB for the molecule 
on the pure metal surface. The electronic configuration is 
close to {dxy)'^ ^id^^)'^ ^{dz2y . The main difference is that 
the partially occupied d^^ orbital close to is now unoc- 
cupied. However, the dj;2_y2 orbital is partially occupied 
such that the ratio of spin-up and spin-down charge den- 
sity is about 2.5 and a clear decision whether the system 
is in a Fe^+ or Fe^+ configuration is not possible. The 
size of the HOMO-LUMO gap is unaffected by the 01 lig- 
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FIG. 12: (Color online) Density of states of Fe d-states and 
CI (gray) in FeP(Cl) on Cu(OOl) (a) and on O/Cu(001) (b) 
with van der Waals interaction. The position of the molecule 
corresponds to the respective ground state configuration, see 
Sec. |V A| Arrows denote the spin-up and spin-down configu- 
ration, respectively. 
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TABLE II: Calculated Fe-surface distances c/pe-surface and 
bond length of Fe and ligand atom (iFe-iigand for FeP on 
missing-row position \f2 x 2v^R45° O/Cu(001) (named 
FeP/O/Cu(001) in the table) and on the hollow site position 
on Cu(OOl). In addition the Fe-N distance dpe-N is given. 
In all calculations van der Waals forces employing Grimme 
correction have been taken into account. 



System 


FeP/Cu(001) FeP/O/Cu(001) 


Ligand 


CI 


CI 


C^Fe — surface (A) 


2.64 


3.08 3.16 


C^Fe-ligand (A) 


2.29 


2.26 1.68 


C^Fe-N (A) 


2.04 


2.02 2.09 



and, but the splitting between the unoccupied d density 
of states becomes larger, being 1.63 eV instead of 1.38 eV 

In the experiment 



without the extra ligand, see Fig. [12 
this would express itself in a broadening of the L3 XAS. 
Furthermore, the pre-edge features are expected to be- 
come smaller, because the majority drj.i_yi orbital forms 
now the LUMO and has the same energy as the minority 



6^7^ orbitals, see discussion in Sec. |VB 

Assuming that an O atom is attached to the molecule 
(similar to the case of CI Fig. 11 ), or more precisely to the 



Fe atom, the scenario is very different. On the one hand 
the structural changes are similar to the previous case, 
i.e., the bonding strength between surface and molecule 
decreases due to the ligand. In the presence of an O 
ligand the molecule-surface distance amounts to 3.16 A, 
whereas the bond length between Fe-atom and ligand is 
only 1.68 A. The small bond length is comparable to the 
findings in bulk FeO, in which the Fe-0 distance is about 
1.62 A."^ On the other hand the oxygen atom influences 
the magnetic properties of the molecule much more than 
the above discussed CI atom. Here, the spin moment 
of the Fe atom is 3.50 /iB, c/. Tab. |III| In contrast to 
the previous case the Fe atom is now in the S = 2-like 
spin-state with /n^ = 4.66/1.20, i.e., only one spin- 
down (mostly character) orbital remains occupied, 
see Fig.[T3j The small bond length between Fe and O 
as well as the tendency of Fe to oxidize lead to strong 
hybridization effects between molecule and ligand, i.e., of 
O p orbitals (gray shaded in Fig.[l3| with Fe d^-^ and d^^ 
orbitals. Other than for CI ligands O strongly influences 
the HOMO-LUMO gap. A p-d hybrid orbital appears 
directly at in the spin-up channel, whereas the large 
gap in the spin-down channel remains nearly unchanged. 
Summarizing, the binding of CI or O ligands to the Fe 
center of the porphyrin molecule stabilizes the high-spin 
state. The occurrence of the high-spin state goes along 
with an increase of the Fe-N distance. In contrast to the 
previously discussed systems without ligands here the Fe- 
N bond length lies between 2.02 to 2.09 A, see Tab.jn) 
These observations agree with findings for Fe porphyrin 
molecules in the gas phase and on graphene'^. 
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FIG. 13: (Color online) Calculated density of states of Fe and 
O (gray) in (FeP(O)) on ^ x 2^R45°O/Cu(001) including 
van der Waals interaction. Arrows denote the spin-up and 
spin-down configuration, respectively. 



A. Comparison between experiment and theory 

1. XAS vs. DOS 

In order to compare the theoretical results to X-ray 
absorption experiments the sequence of the unoccupied 
d states can be related to the features of the XAS at the 
Fe L3 edge obtained from linear polarized light. Different 
parts of the density of states will be visible for different 
incidence angles of the photon beam and polarization 
orientations. For normal incidence d orbitals in the {x — 
y) plane are scanned and in case of grazing incidence 
mostly orbitals with z components contribute, see Fig. [3] 
In the following we focus on the Fe L3 edge only and for 
a discussion the energy levels of the unoccupied d states 



are included in the graph, see Fig. [14 

Although we are aware of the fact that a direct com- 
parison of density of states and XAS is difficult because 
of missing core-hole effects in theory, it gives at least in- 
formation of the character of the XAS peaks. The XAS 
of Fe OEP on Cu(OOl) shows a broad Fe L3 edge signal 
for grazing {0 = 25° (horizontal)) incidence. This corre- 
sponds to the broadened d^^ and orbitals which we ob- 
tained from the present DFT calculations, see Fig. 14 ^a). 
In case of normal incidence two peaks occur at 706.8 and 
708 eV. This is in reasonable agreement with the sepa- 
ration of the two d^2_,,2 orbitals which is 1.33 eV. The 



pre-peak feature in Fig. 14 is more pronounced for nor- 
mal incidence, i.e., the electrons are excited to d orbitals 
which are (partially) in-plane oriented. This corresponds 
to the d^^ states which form the LUMO orbital. Here, 
the overall agreement between the calculation and the 
experiment is good both in view of orbital character and 
energetic distance between the peaks, cf. Fig. [14 



As discussed in Sec.[IV| the X-ray absorption spectra 
turned out to be very different on the oxidized surface, 
see also Fig.[l4[b). Here, the intensity of the Fe L3 edge 
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spectrum is rather independent from the photon inci- 
dence angle 6 and the peaks are sharper compared to 
Fe OEP on the Cu(OOl) surface, which coincides with 
the fact the calculated energy levels are less spread com- 
pared to the FeP molecule on the plain Cu(OOl) surface, 
see Fig.[T4] The d^^ and orbitals are actually located 
between the two dx2_y2 orbitals such that the XAS and 
the calculated density of states are in qualitative agree- 
ment. However, in case of the oxidized surface the dis- 
tance between the two in-plane polarized levels does not 
match perfectly with the distance between the two main 
features of the L3 edge which may be related to core-hole 
effects. 

The situation becomes more complex if the stabiliz- 
ing ligand is CI instead of pyridine, cf. Fig. [151 



OEP(Cl) on Cu(OOl) the intensity for grazing incidence 
angles is now larger and lies basically between the fea- 
tures obtained at normal incidence of the photon beam 
such that the highest intensity is observed for grazing in- 
cidence {0 = 70°). A comparison of the XAS of the Fe 
L3 edge of FeP (CI) /Cu (001) with the energy levels of the 
unoccupied Fe d states shows that the position of the two 
dx2-y2 levels is in reasonable agreement with experiment, 
whereas the angular dependent intensity of the XAS for 
grazing incidence barely fits with the calculated d states 
of FeP/Cu(001), see Fig[T5| Using the level spectra of 
FeP(Cl)/Cu(001) instead, the high intensity of the XAS 
around 708 eV is quite well reproduced by the broad dis- 
tribution of d-j^ and d^^ states. However, the 2nd dx2_y2 
level is close to 706.5 eV and therefore, about 0.5 eV too 
low. These results can be viewed as an indication that 
CI has been only partially removed and the coexistence 
of Fe OEP with and without CI ligand is reflected in the 
XAS of the Fe L3 edge. 

While the XAS of the Fe L3 edge of Fe OEP on Cu(OOl) 
shows a significant dependence on the ligand used to sta- 
bilize the molecule, this is not observed for the Fe OEP 
molecule on the oxidized surface, cf. Fig. [14] and [15] In- 
dependent from the preparation of the Fe OEP molecule 
the angular dependence of the XAS is quite similar. From 
a comparison of the XAS of the Fe L3 edge with the 
calculated energy levels, it is obvious that the spectrum 
is already weh described by FeP on O/Cu(001). The 
broadening of the d^ and orbitals which occurs due 
to CI does not fit with the experimental results. Hence, 
the influence of CI seems to be less pronounced in this 
case. However, more information is needed to underpin 
the findings from the comparison of XAS and calculated 
energy spectra. This is done in the next section by ana- 
lyzing the magnetic properties of the four systems. 



2. Effective moments and magnetic anisotropy 

The measured XMCD of the Fe L3 edge of Fe OEP 
shows a crucial dependence on the photon incidence angle 
^, i.e., the signal is very tiny for normal incidence but 
large for smaller angles. This is an indication that either 
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FIG. 14: (Color online) XAS of Fe L3 edge of Fe OEP on 
Cu(OOl) (a) and v^x 2x/2 R45°O/Cu(001) (b) measured with 
linear polarized light. The calculated energy levels of the un- 
occupied Fe 3d orbitals are shown as horizontal lines. The 
percentage of the d character is given in brackets. The per- 
centage of the d character is given in brackets. Black lines 
mark the XAS with electric field mostly oriented perpendic- 
ular to the surface {6 — 25° horizontal) and the bright line 
denotes the XAS obtained for — 25° (vertical), where the 
electric field vector is parallel to the surface. 



the dipole term l{Tz) is large or that the system has a 
large magneto-crystalline anisotropy. Both quantities are 
not directly accessible in XAS and XMCD experiments. 
The spin magnetic moment which is obtained from the 
sum-rules is an effective spin moment meff and is related 
to the actual spin moment by 



^eff 



(4) 



For bulk systems the dipolar term is usually small and 
can be neglected, however, it may be significantly large 
in hybrid systems consisting of molecules and surfaces.^ 
To determine 7{Tz) from the ah initio calculated density 
matrix we followed the method of van der Laan.^^ This 
model allows the prediction of the effective moment for 
different photon angles assuming the system is magnet- 
ically fully saturated, see Sec.|llJ The calculated 7{Tz) 
term for FeP on Cu(OOl) is given in Fig. 16 by the black 
line and symbols. Indeed the calculated dipolar term is 
large, i.e., of the same order as the spin moment itself. In 
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FIG. 15: (Color online) XAS of Fe L3 edge of OEP(Cl) on 
Cu(OOl) (a) and 2^/2 R45°O/Cu(001) (b) measured with 
linear polarized light. The calculated energy levels of the 
unoccupied Fe 3d orbit als are shown as horizontal lines. The 
percentage of the d character is given in brackets. Black lines 
mark the XAS for normal incidence {0 = 0°) for which the 
electric field lies in the plane and the bright line denotes the 
XAS obtained for grazing incidence {0 = 70°) with the electric 
field being mostly perpendicular to the surface. 



saturation the calculated effective spin moment for graz- 
ing incidence with = 65° {0 = 90°) should be 2.45 //b 
(3.0 /iB)- The small negative value obtained for ^ = 0° 
is related to the limited accuracy of the determination 
of the dipolar term and has no physical meaning. One 
reason may be related to the fact that the d^^ and d^^ 
orbitals, which give the main contribution at ^ = 0°, are 
less sharp than the in-plane oriented orbitals. The ac- 
tual values for meff obtained from X-ray absorption spec- 
troscopy are much smaller namely 1.06 /iB at the magic 
angle (see Tab. Ill), because the maximum magnetic field, 



which was available during the measurement, was 5.9 T 
which was not sufficient to saturate the sample, see Fig.|5| 
The occurrence of large dipolar contributions for systems 
with reduced symmetry has already been predicted by 
Crocombette et a/.*^ who investigated the dipolar term 
for 3d transition metals with different local symmetry. 
For FeP on Cu(OOl) it turned out that the calculated 
values, which reach from -2 to +l/iB, are of the same 
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FIG. 16: (Color online) Effective spin magnetic moments of 
Fe in FeP/Cu(001) depending on the incidence angle of the 
photon beam d. The theoretical effective moments are ob- 
tained from the calculated density matrix. The black line 
denotes the calculated "l^z) term and the dashed lines marks 
the magic angle where in saturation vfis — rrieff. 



order as the spin moment = 1.95 /iB- This agrees 
qualitatively with recent findings by Stepanow et al. for 
Cu phthalocyanine on Ag(OOl).— 

In the presence of oxygen on the Cu(OOl) surface the 
dipolar term, or more precisely, its angular dependence 
is almost identical (not shown here). The magnetic spin 
moment for = 90° amounts to 3 /iB as in the case with- 
out O. This behavior is expected because the occupa- 
tion of the Fe d levels is very similar in both cases, see 
Sec. |VB[ which is reflected in the density matrices. In 
this case the experimental results for Fe OEP (Py) on 
the oxidized surface agree well with the theoretically ob- 



TABLE III: Effective spin moments of Fe obtained from sum- 
rule analysis and calculated spin moments of Fe (includ- 
ing van der Waals forces) in Fe porphyrin on O/Cu(001) 
and Cu(OOl) together with the calculated ratio of spin-up 
and spin-down charge of the Fe d orbitals and the magnetic 
anisotropy energy (MAE). Magnetic moments are given in /jb 
and energies in meV. The experimentally determined moment 
are given at the magic angle, whereby not all systems were in 
saturation, see text. 



System 




^Fe/^Fe 




MAE 




(/is) 


(%) 




(meV) 


Experiment: 










FeOEP(Cl)/Cu 


>1.01±0.15 


13 






FeOEP(Cl)/0/Cu 


>2.32±0.30 


8 






FeOEP/Cu 


>1.06±0.15 


26 






FeOEP/O/Cu 


2.24zb0.30 


12 






Theory: 










FeP/Cu 


1.93 




3.96/2.07 


0.11 


FeP/O/Cu 


1.99 




3.99/2.05 


-0.10 


FeP(Cl)/Cu 


2.50 




4.18/1.69 


-0.09 


FeP(Cl)/0/Cu 


2.62 




4.21/1.64 


-0.10 


FeP(0)/0/Cu 


3.50 




4.66/1.20 


0.45 
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tained effective spin moments, see Fig.[T6j because the 
sample is almost saturated for the maximal applied field 
of 5.9 T. At the magic angle both experiment and the- 
ory give a spin magnetic moment of 2 /iB • The good 
agreement between experimental and calculated effective 
spin moments at the magic angle suggests the assump- 
tion that the angle dependent changes in the XMCD are 
not related to a large MAE. 

The dependence on the photon incidence angle is quite 
different for the Fe OEP molecules which have been sta- 
bilized with CI ligands. The effective moment obtained 
from experiment at the magic angle is 2.70 /iB whereby 
corresponding to F ig. [5] the sample is not fully saturated, 
see squares in Fig.[T7| and the calculated meff without 
ligand clearly underestimates the magnetic moment, cf. 



dashed line in Fig.[T7| However, the calculations for FeP 
with CI result in an effective spin moment which is larger 
than the experimental value whereas without CI the effec- 
tive moment is below the experiment. This underpins our 
previous observations that CI has not been completely 
removed, cf. Sec. |III| The remaining CI has obviously a 
significant infiuence on the magnetic configuration. Oxy- 
gen ligands lead to a completely different angular depen- 
dence and can be ruled out as possible cause for the dif- 
ferent magnetic structures of Py- and Cl-ligated Fe OEP 
molecules on Cu(OOl) and O/Cu(001). 

However, a huge dipolar term is only one possi- 
ble reason for the strong angular dependence of the 
XMCD (Fig.[6]), the other possible explanation would 
be a large magnetic anisotropy of the system. Here, 
we have used the method proposed by Wang^^ to ob- 
tain the anisotropy energy from 2nd order perturba- 
tion theory. The calculated magnetic anisotropy ener- 
gies obtained for the porphyrin molecule on the plain 
and the oxidized Cu(OOl) surface are -0.096 (0.114) eV 
with (without) O (minus sign means out-of-plane easy 
axis). The MAEs are by a factor of 10 smaller than 
the values obtained by Wang et at. for Fe phthalocya- 
nine molecules, which have a structure very similar to 
the FeP molecules investigated here. Wang et al. calcu- 
lated for the gas phase molecule an in-plane anisotropy of 
1.18 meV.^^ However, these calculations were performed 
within the GGA and therefore, the distance between oc- 
cupied and unoccupied states was underestimated, which 
in turn leads to larger anisotropy values, see Sec.jlTj The 
MAE obtained for the FeP molecules on Cu(OOl) is also 
smaller than for free Jahn- Teller-distorted Feis clusters, 
in which the calculated magneto-crystalline anisotropy 
is about 0.6meV/Fe atom.^^ Also in Fe-Pt nanoparticles 
which have been capped with Cu the magneto-crystalline 
anisotropy amounts only to 0.5 meV/atom.!^ Compared 
to that the calculated MAE for FeP on the Cu(OOl) sur- 
face are quite small and the strong changes in the XMCD 
spectra cannot be explained from the magneto-crystalline 
anisotropy, but may be related to the dipolar term and 
its angular dependence. Unfortunately, due to the large 
super cells the accuracy of the calculations is limited such 
that the determination of MAEs is already at the edge of 




FIG. 17: (Color online) Calculated effective magnetic spin 
moments of Fe and dipolar term 7{Tz) depending on the inci- 
dence angle of the photon beam in the presence of different 
ligands at the Fe-atom in FeP/V2 x 2^/2 R457O/Cu(001). 
Diamonds denote an 0-ligand and circles correspond to a CI 
ligand. The dashed curve corresponds to the calculated meff 
without ligands. The green squares mark the experimental 
data for Fe OEP(Cl)O/Cu(001). 



the accuracy. This means that for small anisotropics the 
determination of the easy axis is difficult. Except for FeP 
on Cu(OOl) the calculations indicate an out-of-plane easy 
axis, which is not supported by the experimental findings, 
see Fig.|5j It should be mentioned that the infiuence of 
the van der Waals forces on the magnetic anisotropy en- 
ergy is negligible in our case, since for the present system 
the magnetic configuration remains nearly unchanged if 
van der Waals interaction is taken into account. Anal- 
ogous to the previous investigations without ligands, we 
calculated the effective spin moment and the dipolar term 
{Tz). In the presence of CI ligands, the dipolar term 
would be even more important than before. Our cal- 
culations of Tz reveal that meff strongly depends on ^, 
i.e., the dipolar term is huge being -2.6 /iB for normal 
incidence and 1.4/iB at 90°, cf. Fig.[T7l The magnetic 
anisotropy obtained from 2nd order perturbation theory 
for the Cl-ligated FeP on the O/Cu(001) (Cu(OOl)) sur- 
face is -0.10 eV (-0.09 eV), i.e., it is of the same order of 
magnitude as for the systems without CI. 

Although the O ligand case is somewhat artificial, we 
have also determined the dipolar term and the effective 
magnetic moment. The angular dependence of Tz ranges 
now from -0.63 to 1.25 /ie, which is small compared to 
the previous cases without ligand and Cl-ligated FeP 
molecules, see Fig.[T7[ However, the most striking differ- 
ence is the sign of the Tz term. The average dipolar term 
is positive being 1.26 /iB and has the same sign as in the 
magnetic spin- moment. The difference arises from differ- 
ent occupation of the spin-down orbitals. Here, mainly 
dz'2 and d^^ orbitals are occupied, whereas in case of CI 
the minority channel has mainly dxy character. In this 
case the magnetic anisotropy from 2nd order perturba- 
tion theory amounts to 0.45 meV which is by a factor of 
4 larger than for all other cases. 
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VII. CONCLUSION 

We have presented a combined experimental and theo- 
retical study of the magnetic and electronic properties of 
Fe GEP molecules on nonmagnetic substrates. Fe GEP 
molecules on metallic Cu(OOl) and an oxidized Cu(OOl) 
surface have been investigated in view of the spin-state, 
the MAE, and the influence of ligands on the magnetic 
and electronic properties of the molecule-surface hybrid 
system. XAS and XMCD have been measured for differ- 
ent incidence angles of the photon beam. In atmosphere 
Fe GEP molecules are stabilized by ligands which bond 
to the Fe atom. Here we have used two different samples 
one with atomic CI and one with pyridine, which turned 
out to have a significant impact on the XAS. 

The spin state and the adsorption position of 
the molecules on fee Cu(OOl) and on the a/2 x 
2^/2 R45°G/Cu(001) surface have been obtained from 
DFT calculations. To account for correlation effects, van 
der Waals forces were included in the approximation of 
Grimme. Here, the inclusion of the long-range interac- 
tions basically reduces the molecule-surface distance by 
0.5 A. In contrast to adsorption on magnetic surfaces^ 
the molecules are weakly bonded, i.e., the hybridization 
with the surface is small, however, density of states as 
well as XAS show distinct differences, which are related 
to the structure of the substrate. Also the interaction 
mechanism with the surface changes. In case of Cu(OOl) 
the hybridization effects and transfer effects indicate that 
interaction between Fe atom and surface is indirect via 
the N atoms, which is usually observed on magnetic metal 
surfaces. In the presence of the oxygen layer a significant 
part of the hybridization and interaction with the sur- 
face is mediated by the oxygen atoms. For a quantitative 
comparison of the effective spin moments obtained from 
sum-rule analysis of XAS and calculated spin moments, 
we determined the dipole term which connects the two 
quantities. 

From DFT calculations an intermediate spin state 
with 2 fiB is obtained for FeP molecules on both sub- 



strates, whereas the absorption sites are different, namely 
missing-row position with oxygen and hollow site position 
on the unreconstructed metal surface. The magnetic mo- 
ment and the spin state are in agreement with the exper- 
imental results with pyridine as stabilizing ligand in air, 
whereby only in the case of the oxidized surface mag- 
netic saturation could be approached. The calculated 
magneto-crystalline anisotropy energy is about 0.1 meV 
such that the strong angular dependence of the XMCD 
has to be mainly related to a huge dipolar term, which 
is supported by our calculations. 

If CI has been used to stabilize the Fe GEP molecule, 
XAS and XMCD show a different dependence on the pho- 
ton incidence angle, whereby the effect is more expressed 
on the Cu(OOl) surface. A comparison of the Fe L3 edge 
XAS with the energy levels of the unoccupied Fe 3d or- 
bitals gave rise to the assumption that CI has not been 
completely removed. This was underpinned by the anal- 
ysis of the Tz term and the calculation of the effective 
spin moments. 

Summarizing, the combination of X-ray absorption 
spectroscopy with DFT calculations and, in particular, 
the analysis of the term give an accurate picture of 
the electronic and magnetic structure of Fe porphyrin 
molecules on the nonmagnetic substrate and help to dis- 
tinguish effects from the surface and remaining ligands 
attached to the molecule. 
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